Collective coherent population trapping in a thermal field 
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We analyzed the efficiency of coherent population trapping (CPT) in a superposition of the ground 
states of three-level atoms under the influence of the decoherence process induced by a broadband 
thermal field. We showed that in a single atom there is no perfect CPT when the atomic tran- 
sitions are affected by the thermal field. The perfect CPT may occur when only one of the two 
atomic transitions is affected by the thermal field. In the case when both atomic transitions are 
affected by the thermal field, we demonstrated that regardless of the intensity of the thermal field 
the destructive effect on the CPT can be circumvented by the collective behavior of the atoms. 
An analytic expression was obtained for the populations of the upper atomic levels which can be 
considered as a measure of the level of thermal decoherence. The results show that the collective 
interaction between the atoms can significantly enhance the population trapping in that the popula- 
tion of the upper state decreases with increased number of atoms. The physical origin of this feature 
was explained by the semiclassical dressed atom model of the system. We introduced the concept 
of multiatom collective coherent population trapping by demonstrating the existence of collective 
(entangled) states whose storage capacity is larger than that of the equivalent states of independent 
atoms. 
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I. INTRODUCTION 

The study of atomic coherence effects in multilevel 
atoms is one of the most active area in atomic spec- 
troscopy T, "5, "s"! . Especially, the theory of coherent pop- 
ulation trapping (CPT) in a three-level A-type atom has 
been extensively studied and the phenomenon has been 
observed experimentally in a sodium vapor 0,13, pho- 
toassociation systems BEC (3 and solids [g. The 
CPT results from the formation of a coherent superposi- 
tion of the ground atomic states that is decoupled from 
the external fields and hence referred to as a dark state. 
The particular interest of this phenomenon consists of 
the possibility of storage and coherent manipulation of 
the population in a coherent superposition of the ground 
states of the atoms . These phenomena have re- 

ceived greatly increased experimental attention in recent 
years and experimental techniques have been developed 
which allow a reversible transfer of quantum informa- 
tion from light to the dark state of the atoms jllj . The 
coherent population trapping has also been investigated 
in the context of lasing without inversion ) subrecoil 
laser cooling and a search for materials that display 
a high index of refraction accompanied by vanishing ab- 
sorption [miiiii. 

The atomic coherence effects are sensitive to decoherence. 
In the CPT effect, one source of decoherence is fluctua- 
tions of the laser fields used to create the coherent super- 
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position of the atomic ground states 0| . The fluctua- 
tions redistribute the population among the atomic states 
including the excited atomic states from which it can be 
spontaneously emitted resulting in optical losses. Recent 
investigations of decoherence processes in atomic systems 
have demonstrated that CPT and quantum storage in an 
ensemble of noninteracting atoms are limited primarily 
by different decoherence processes such as atomic colli- 
sions, atom loss and motion of atoms 0,^3- The results 
show an interesting property that in the limit of the total 
number of excitations much smaller than the number of 
atoms, the decoherence rate of the multiatom system is 
of the same order of magnitude as in the single atom, i.e. 
is independent of the number of atoms in the sample. In 
an earlier study, Jyotsna and Agarwal 19] showed that 
the CPT effect in a dense atomic medium is unaffected 
by local-field effects. 

It is well known that the dominant contribution to the 
decoherence processes in the interaction of atoms with 
the electromagnetic field stems from the thermal fluc- 
tuations. They are present in a non-zero temperature 
reservoir to which the atoms are coupled. The fluctu- 
ations cause a pumping of the population stored in the 
dark state into the excited states of the atoms from which 
it can be spontaneously emitted resulting in an increase 
in decoherence. The magnitude of thermal fluctuations 
depends on temperature of the reservoir and determines 
the minimum level of thermal decoherence. 

In this paper we propose a method to suppress the deco- 
herences that occur due to the thermal fluctuations of the 
environmental electromagnetic reservoir at temperature 
T. Essentially, we examine the CPT effect in three-level 
A systems by addressing a practical question: How can 
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one increase the efRciency of trapping and storage of the 
population in the presence of thermal decoherence. In 
particular, we will investigate limits to the efficiency of 
the CPT effect in a single atom and next will explore 
the role of multiatom collective behavior in the reduc- 
tion of the single-atom decoherence rate induced by the 
thermal field. The dipole-dipole interactions between the 
atoms will not be taken into account here assuming lower 
atomic densities, so that the collective behavior we con- 
sider stems entirely from the mutual coupling of all the 
atoms with the common radiation field 20J. Employ- 
ing the analytic solution for the density operator of the 
system, we find that in general the single-atom coherent 
population trapping effect, reduced by thermal fluctua- 
tions, can be significantly improved or even completely 
restored when the atoms interact collectively with the 
thermal modes of the reservoir. We are particularly in- 
terested in the manner in which multiatom effects can 
lead to a suppression of thermal decoherence. With ap- 
propriate selection of atomic parameters, we will find 
cases of almost perfect coherent population trapping in 
the presence of the thermal decoherence. Our physical 
interpretation of the results is based on the semiclassi- 
cal dressed atom model of the collective atomic system. 
The collective dressed states of the system are identified, 
and the effect of suppression of the thermal decoherence 
is explained in terms of the increased capacity of these 
states. This is shown to arise from correlation-enhanced 
transition rates among the multiatom dressed states, in 
particular those entering the trapped state. Hence, the 
effects of decoherence by thermal fields may by reverted 
more rapidly. 

II. APPROACH 

The system we consider is an ensemble of A'' identical 
three-level A- type atoms each with excited state |1) and 
two nondegenerate ground states |2) and |3). The atoms 
are driven by two single-mode cw laser fields of Rabi 
frequencies 21^2 and 2Q3 and angular frequencies ujl2 
and significantly different from each other, so that 
each laser is coupled only to one of the allowed transi- 
tions, as shown in Fig. ^ The transitions are associated 
with nonzero dipole moments ^12 and /ii3, and the laser 
fields are detuned from the atomic transition frequen- 
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FIG. 1; Energy-level diagram of a three- level A- type atom 
driven by two laser fields of Rabi frequencies 2i}2 and 20,^. 



cies, such that there is a nonzero two-photon detuning 
A = (wi3— a;i2-l-Wi2— ^L3)/2. The transition |2) ^ |3) is 
forbidden in the electric dipole approximation (/I23 = 0) . 
The exited atoms may decay spontaneously due to the 
zero point fluctuations of the electromagnetic field from 
the state |1) to both ground states |2) and |3) with the 
decay rates 272 and 273, respectively. We assume that 
the atoms are contained in a volume with linear dimen- 
sions that are small compared with the radiation wave- 
lengths, the Dicke model Thus, all atoms experience 
the same Rabi frequencies of the driving fields including 
their phases, and propagation effects are negligible due 
to the small size of the sample. In addition, we assume 
that the atomic transitions are driven by a thermal field 
of the mean photon numbers n2 and 77.3 at the atomic 
transition frequencies ^12 and 0^13, respectively. 
The system is described by the reduced density opera- 
tor, which in the interaction picture and under the usual 
Born-Markov and rotating- wave approximations satisfies 
the master equation 

^ ^ -J^[Hq,p\+"/2C2P + J3^3P, (1) 

where 

Ho = hA{S22 - S33) + h + 

ae{2,3} 

C2P = (1 + ri2)[5'2ip, S'12] +n2[5'i2p, S'21] +H.C., 

C3P = (l + n3)[^3lP,^13]+n3[5l3P,53l]+H.C.. (2) 

Here C2P and C3P are operators representing the damp- 
ing of the atoms via spontaneous emission and Hq is the 
Hamiltonian describing the coupling of the atoms to the 
laser fields. The operators Sa0 are the collective atomic 
operators 

N N 

^«/3 = E^S «:/3 = l,2,3, (3) 

which obey the usual commutation relations 

The master equation allows to obtain equations of 
motion for the expectation value of an arbitrary com- 
bination of the atomic operators. The calculations can 
be performed without much troubles for the simple case 
of a single atom {N — 1) and arbitrary A. However, 
for > 1 the calculation of the expectation value is 
not an easy task. In even the simplest cases of small 
numbers of atoms, the calculations are prohibitively dif- 
ficult due to the enormity of the number of coupled equa- 
tion of motion. Fortunately, for the A = case and 
high field strengths, flk 3> Njk, an approximation tech- 
nique has been developed, which greatly simplifies the 
master equation and thus ables to perform analyti- 
cal calculations of the expectation value of an arbitrary 
combination of the atomic operators. The restriction to 
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the A = case stems from the difHculty in obtaining a 
closed set of equations when the two-photon detuning is 
present 123 A full discussion of the technique is given 
in Refs. [H H HI El IH IH El 113 ■ In the interest 
of brevity only the key results will be given here. The 
technique is implemented by introducing dressed states 
of a single atom, which are obtained by a diagonalization 
of the single-atom interaction Hamiltonian 
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(5) 



The single-atom dressed states are of the form 



l*i>. 

l*3>, 



r!3|2)j 
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71' 



V2n 
1 

V2n 



f^3|3),), 
^^3|3),), 



(6) 



where 57 = is the generalized Rabi frequency. 

The idea of the approximate technique is now to replace 
the collective operators Sap by the collective dressed- 
atom operators 



N 



N 



^a;3 = E^S = El*")^-j(*/3|, a,/3 = 1,2,3, (7) 



and then substitute for Sa0 into the damping terms of the 
master equation Next, we make the unitary trans- 
formation of the density operator 



exp 



■^Hot) pexp ( --^Hot 



where 



Ho = hQ {R22 - R33) ^ hQR, 



(8) 



(9) 



and on carrying out this procedure it is found that cer- 
tain terms in the transformed master equation are slowly 
varying while the others are rapidly oscillating at fre- 
quencies and 2f2. The approximation then consists 
of dropping these rapidly oscillating terms. The master 
equation |^ in the dressed state basis reduces to 



dp 
dt 



where 
To = 

Ti = 

r2 = 



-in[R,, p] + {Toi[R,p, R,] + [i?32P,i?23] 

[R23P, R32]) + ri([i?i2P, i?2i] + [Ri3p, R31]) 

r2([i?2lP, R12] + [R31P, R13]) + H.C.}, (10) 



- 73(1 + 2^3) [ 



^3 



{72(l + n2)[r!3/r!]- 



V2n 

73(l-|-n3)[r!2/f^]'}, 



{j2n2[^3/^? +I3n3p2/n]^}, 



(11) 



are the transition rates between the single-atom dressed 
states. 

Using the approximate master equation, it is straightfor- 
ward to obtain a simple analytical solution for the steady- 
state density operator of the system. The solution can 
be written in the form 



Ps ^ Z ^cxp[-^i?ii]. 



where 



^ = ln 



fl^n2 + ?7ri2?T-3 



f7^(l + n2)+ 77^^2(1 + "•3) 



(12) 



(13) 



and rj = 73/72. The parameter Z is the normalization 
constant such that Tr{/5s} = 1. It is easily verified that ^ 
is always negative independent of the parameters used 
and approaches zero when h2 and/or 7x3 go to infinity. 
The solution (O was obtained in Refs. pi I2I Isq. 
and some applications are discussed there in details. In 
Ref. [13 , the solution has been used to investigate differ- 
ent control schemes for collective systems of three-level 
atoms. In this paper, we focus on the competition be- 
tween thermal fiuctuations and the collective effects that 
can lead to collective population trapping. 
The steady-state solution (|12f) enables to calculate any 
statistical moment of the diagonal elements Raa, and 
thus population distributions between atomic states. In 
particular, an fc-th order moment of Rn (expectation 
value of a product of k operators i?ii), is of the form 



{R^,,), = {-lfZ-^^Z, fc=l,2,..., (14) 
and the first order statistical moments of i?22 and i?33 are 



(-R22). = {R33)s = [N-{Rn)s]/2, 



where 



Z 



A^ + 2-(jV + l)e^-e-^(^+i) 
(l-eC)(l-c-«) ■ 



(15) 



(16) 



One can easily show from (|12|l that the steady-state off- 
diagonal elements Ra^ {a ^ (3) equal zero. Note from 
Eq. (dS)) that aU the non-zero expectation values can be 
represented in terms of {Rii)s- The steady-state solu- 
tions are to be used in the forthcoming treatment of the 
coherent population trapping in a multiatom system. 



III. COHERENT POPULATION TRAPPING 

Before we proceed to the detailed analysis of the mul- 
tiatom trapping effect, we briefly investigate the trap- 
ping behavior of single atoms in the presence of thermal 
fluctuations. In this way we may see what restrictions 
are brought by the thermal fluctuations for the trapping 
phenomenon and how they are related to the coherent 
driving process. Coherent population trapping effect in 
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FIG. 2: Stationary population of the upper state as 
a function of the two-photon detuning A for 72 = 73 = 7, 
f22 = ^3 = 57 and different n: n = (sohd hue), n = 0.5 
(dashed hue), n = 2 (short dashed hue). 



a system of three-level atoms may be monitored experi- 
mentally in terms of the intensity of the fluorescence light 
emitted 0, IE IE S 13 ■ It is manifested by the disappear- 
ance of the fluorescence which, on the other side, mani- 
fests the vanishing of the population of the upper atomic 
states Therefore, we will consider first the effect of 
the thermal field on the so-called transparency window, 
i.e. the dependence of the stationary population on 
the two-photon detuning. Next, using the stationary so- 
lution H12|) , we will find the analytical expression for the 
population at the two-photon resonance, A = 0, and will 
analyze how one could reduce the destructive effect of 
the thermal field on the minimum of the population at 
A = 0. 

Fig. ^ illustrates the stationary population as a 
function of the two- photon detuning A. We have ob- 
tained the population by solving numerically the master 
equation ^ for iV = 1. It is seen that in the absence of 
the thermal field, n = 0, there is perfect CPT observed 
at A = 0. When the atom is in the thermal field equally 
affecting both transitions, the CPT effect is reduced and 
the thermal field washes out the transparency window as 
n ^ 1. Thus, the thermal field has a destructive effect 
on the CPT, because the thermal field is an incoherent 
field with random fluctuations that destroy the coherent 
process induced by the laser fields. 

The variation with n of the minimum of the upper state 
population at A — can be analyzed explicitly using 
Eq. © which for = 1 and together with the steady- 
state solution (|12|l gives a simple analytical expression 
for pf^ ill the form 

Pn-^((i?22). + (i?33>s) = (17) 



First, we note from Eq. (|17() that the population dis- 
tribution between the atomic states is determined solely 
by the parameter ^. Clearly, the population distribution 
and consequently the trapping effect depend on several 
parameters such as the laser intensity, spontaneous emis- 
sion rates, and mean number of thermal photons. We 
can call the parameter ^ as a measure of efficiency of the 
CPT effect. 



Here the efficiency of the CPT is examined in various in- 
tensity regimes of the coherent fields for equal and also 
unequal average numbers of thermal photons. The aver- 
age numbers can be made unequal by a suitable choice 
of bandwidth of the thermal field. The selective excita- 
tion of the atomic transitions can be realized in practice 
by applying a finite bandwidth multimode thermal field 
whose bandwidth is much smaller than the splitting of 
the lower atomic levels, but large compared to the nat- 
ural linewidths of the atomic transitions to satisfy the 
Markov approximation used in the derivation of the mas- 
ter equation. 

In the limit of 73.2 = that the thermal fluctuations affect 
only the |1) |3) transition, the parameter ^ reduces to 



(18) 



The parameter ^ does not change substantially with the 
Rabi frequencies unless fl^ is much larger than the Rabi 
frequency Q2 of the other transition. In the very strong- 
field regime of fig ^ V^2i^ + "-3); the parameter ^ ap- 
proaches the limit of ^ —> —00. This minimum value is 
that one which leads to vanishing of the population of the 
upper atomic state, because lim^_,_oo Pii = 0. This pre- 
dicts that perfect coherent population trapping can be 
observed even in the presence of thermal decoherence on 
one of the two atomic transitions, which is in contrast to 
the result of |^. However, it requires that the transition 
influenced by the decoherence is simultaneously driven by 
a strong laser fleld. It can be understood rather easily. 
For a large Rabi frequency fig, the coherent processes 
on the |1) |3) transition dominate over the incoherent 
thermal processes resulting in perfect transparency. 
Various other intensity regimes can also be distinguished. 
If 7^ the parameter ^ can depend entirely on n2 
or 77,3 depending on the ratio fl3/fl2. For instance, when 
n2^^ S> r/naf^l, we flnd that 



e = ln 



n2 



l + n2 



(19) 



This predicts that the coherent population trapping de- 
pends entirely on the thermal fluctuations at the weakly 
driven |1) — s- |2) transition. In the opposite limit of 
7777,3512 3> 77,2^13, the parameter ^ now depends entirely 
on 77,3. Thus, the driving fields are relatively efficient in 
controlling decoherence in a single atom. Again, it can 
be interpreted as caused by coherent processes that dom- 
inate incoherent thermal processes on the strongly driven 
transition. This also shows that the suppression of the 
thermal decoherence in a single atom is limited to the 
level set by the lowest thermal fluctuations affecting the 
atomic transitions. 

In the case when the thermal fleld equally contributes to 
both atomic transitions, 77,2 = ^^3 = we have 



^ = In 



(20) 
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FIG. 3: Single atom dressed states and possible transitions 
with the rates ro,ri and 



case the transparency window shows a strong sensitivity 
to the Rabi frequencies and is observed only in the limit 
of very weak driving fields. 

The limit set in single atoms by temperature of the reser- 
voir can be circumvented to improve the efficiency of the 
CPT effect if one considers multiatom collective systems 
in which interatomic interactions can create collective 
states of a significantly enhanced storage capacity com- 
pared with the capacity of the corresponding states of 
individual atoms. 



IV. COLLECTIVE TRAPPING STATES 



independent of the Rabi frequencies and the spontaneous 
emission rates. Obviously, the trapping effect is reduced 
regardless of how strong are the Rabi frequencies of the 
laser fields relative to the thermal fluctuations. In other 
words, there is no possibility of obtaining perfect popu- 
lation trapping or a control of the decoherence level in 
a single atom when both transitions are equally affected 
by the thermal field. A qualitative understanding of this 
effect can be obtained in terms of the transition rates 
Figure 13 shows the single-atom dressed states and 
the transition rates F. One can see from the figure that 
the population flows into the state with the rate 

Fi, and is removed from this state with the rate F2. The 
state is a linear superposition of only the ground 

states of the atom that it is the trapping (dark) state. 
Therefore, we can call the rate F2 a decoherence rate, as 
it transfers the population from the dark state to the up- 
per state |1) from which it can be spontaneously radiated 
resulting in an increase in decoherence and optical losses. 
Only in the absence of the thermal field, n2 = = 
the transition rate F2 = 0. Evidently, the CPT effect de- 
pends crucially on F2 , and therefore the key to maintain 
a large efficiency of the CPT is to make F2 as small as 
possible. It can be done when the thermal field unequally 
affects the atomic transitions, i.e. when the number of 
thermal photons affecting one of the transitions is dif- 
ferent than on the other transition. For example, when 
n2 <C na, the rate F2 can be made small, proportional 
to n2, by changing the ratio ^■i/^l2. It is easily to see 
from 1)11(1 that in the case of ^ £72, the rate F2 is only 
of the order of n2 despite the fact that there is a large 
number of thermal photons present on the |1) — 13) transi- 
tion. When the thermal field equally contributes to both 
atomic transitions, ^2^2 = 73"-3 = 7^1, and from Eq. Hll|) 
we find that F2 = ^7/2 independent of the Rabi frequen- 
cies of the laser fields. This is the smallest decoherence 
rate one can achieve in the single atom interacting with a 
thermal field that equally affects the atomic transitions. 
The limit is set by the number of photons n that, on the 
other hand, depends on temperature of the reservoir. An 
improvement of the CPT effect in the A— type system 
with asymmetric spontaneous decay rates has been pre- 
dicted in the absence of the thermal field i but in this 



The effects described in Section III can be seen in di- 
lute atomic gases where the interatomic interactions are 
not important. However, a more interesting situation 
emerges as we have considered here atomic samples where 
radiative interactions between the atoms can lead to a 
collective (entangled) behavior of the atoms. Here, we 
include the multiatom effects and calculate the popula- 
tion as a function of the number of atoms and the 
number of thermal photons. 

The upper state population p\i can be evaluated using 
Eq. lO which, together with the steady-state solution 
(|12ll . gives the analytical expression for in terms of ^ 
and N as 



(l-e«)2 J- ^ ' 

In Fig. ^ we present a three-dimensional plot which 
shows that in the absence of the thermal field, i.e. n = 0, 
the stationary population pl^ is equal to zero indepen- 
dent of the number of atoms. Thus, for fi — O the collec- 
tive behavior of the atoms does not affect the trapping 
effect. The presence of the thermal field has a destruc- 
tive effect on the trapping phenomenon that the popu- 




FIG. 4: The upper-state population pli/N as a function of 
n and N for n2 = = n and moderate numbers of atoms. 
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FIG. 5: The upper-state population pli/N as a function of 
n2 ~ — n for different numbers of atoms: A'^ = 10 (solid 
line), N = 100 (dashed line), N = 1000 (short dashed line). 



lation in the upper state is no longer zero and increases 
with increasing number of thermal photons n. However, 
the rate of the increase of the population decreases with 
increasing number of atoms N such that for a suitably 
large N the population pli may remain very small even 
for large n. In other words, the thermal decoherence 
decreases with increasing number of atoms. Thus, the 
collective interactions are relatively efficient in suppres- 
sion of thermal decoherence such that the atoms may 
remain in their ground states even in the presence of the 
thermal decoherence. This is a suprising result as one 
might expect that decoherence should increase with the 
increasing number of atoms. 

Figure shows the population pl^ as a function of h2 — 
= n for different numbers of atoms. Here we see that 
the rate of the increase of the population decreases with 



J 



A''. For a small number of atoms, the population satu- 
rates quickly with n. But, for a large number of atoms, 
a much stronger thermal field is required to reach satu- 
ration. In other words, the collective population stored 
in the ground states is less affected by the thermal fluc- 
tuations than for the case of independent atoms. As a 
consequence, one has a practical scheme to reduce ther- 
mal decoherence and preserve CPT in the thermal field. 



In order to obtain an insight into the physical origin of 
the reduction of thermal decoherence and the improve- 
ment of the population trapping, we examine the energy 
structure of the collective system. In general, in the 
absence of the driving fields, the system can be repre- 
sented in terms of collective symmetric and antisymmet- 
ric states. However, in the case of N identical atoms 
contained in a volume with linear dimensions that are 
small compared with the radiation wavelengths, only the 
symmetric states couples to external driving fields. The 
dipole-dipole interactions between the atoms lead to a 
shift of these states from the laser resonance Thus, 
here the Rabi frequencies should be larger than the shift 
caused by the dipole-dipole effects, i.e. the latter can be 
neglected. The antisymmetric states do not participate 
in the dynamics of the small sample system |21| . There- 
fore, we may limit the dynamics to only those involving 
the symmetric states. Moreover, only the lowest in en- 
ergy symmetric states are of interest in the analysis of the 
collective population trapping. We therefore consider the 
lowest energy states defined as 



|3) 
|2> 



|12) 
etc.. 



^ |22) . 
= |222) = 



|3i,32, . . . ,3jv), 



-4 N 



|3i, . . . , 2i, . . . , 3n), 

1=1 

N 

|3i, . . . , li, . . . , 3n), 

1=1 



1=1 

-i N 



/ , |3l, . . . , 2i, . . . , 2j, . . . , 3Ar), 
i<j=l 
N 

|3i, . . . , li, . . . , 2j, . . . , 3Ar), 

/ , |3i, . . . , 2i, . . . , 2j, . . . , 2fe 

i<j<fc=l 



(22) 



where the binomial coefficients are the normalization con- 



stants. The states (I22f) are superpositions of single- atom 
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FIG. 6: Energy-level structure of noninteracting collective 
states of the A*' atom system. 

product states \m)i (S) \n)j (g) . . . (g) \k)i that are symmetric 
under the exchange of any pair of atoms. For example, 
the state |2) is a linear superposition of the product states 
in which atom i is in the state |2),; and the remaining 

— 1 atoms are in their states \3)j. 
If we now allow the atoms to interact with the laser fields, 
each state |2'^) couples to the first excited states |12*''~^) 



and 1 2*^) with the Rabi frequencies and f^s, respec- 
tively. Figurel^shows the collective symmetric states and 
possible couplings of the two laser fields. As we have al- 
ready mentioned, we limit the presentation to the lowest 
energy levels which will be mixed together by the inter- 
action leading to a ground collective dressed state, which 
is of the main interest here. The lowest energy state is 
the product state |3) = |3i, 82, . . . , Sat). Each succeeding 
state 12*^) is of energy higher by successive increments of 
S — UJ13 — uJi2- Similarly, each succeeding state 112*^) is 
of energy higher by successive increments of (S. It is in- 
teresting to note that the rotating-wave approximation, 
which we are assuming to be valid, ignores coupling of 
states which differ in excitation by two and higher. In 
other words, the laser fields couple only the neighboring 
ground states through the first excited states. It forms 
a two-dimensional chain of A configurations. With the 
state ordering |3), |1), |2), |12), |22), . . . , |2^), correspond- 
ing to the path of successive excitations of the states 12*^) 
by the laser fields, the interaction Hamiltonian Hq can 
be expressed as an infinite tridiagonal matrix 



Ho/h : 



-iVA 







-(A^- 1)A 








-A^A 




n^y/N - 1 

-(A^- 1)A 







^22^2 
-A^A 



(23) 



It is interesting to note that the matrix element of the 
|3) — |1) transition coupled by the Rabi frequency fl^ 
is enhanced by a factor \/N and the magnitude of the 
matrix elements of the successive transitions coupled by 
the same field decreases along the path to the state |2^). 
On the other hand, the matrix element of the |2) — |1) 
transition coupled by the Rabi frequency Q2 is the same 
as in the single atom case, but the the magnitude of the 
matrix elements of the successive transitions coupled by 
the same field increases as Vk when one moves along the 
excitation path to the state |2^). Thus, the coupling 
strength of the lasers to the atoms is transferred from 
one field to the other as one moves along the path of 
excitations from |3) to |2^). 

We now proceed to diagonalize the matrix H23|) which will 
result in collective dressed states. The diagonalization 
is performed by solving Schrodinger's time- independent 
equation in the form 

{Ho - Xnl) \Dn)N = 0, (24) 

where / is the identity matrix and \Dn)N is an eigenvec- 
tor. 



Substituting Eq. (|23(l into Eq. H24|l yields the eigenvalue 
equation 

(A„ -I- NA){[Xr, + {N- 1) A] (•••) + njiXr, 
+ iN - 1)A](. • • )} + ^ImiXn + NA)i- ■ ■ ) 
+rtl{N-l)i---)] = (25) 

where the (• • • ) refers to terms of which the explicit form 
is not needed apart from that those are polynomial func- 
tions of A„. It is easily to show that in the case of A = 0, 
the eigenvalue equation reduces to 

A„ [A„ (• • • ) + ^IXn (•••) + f^3^ (■••)]= (26) 

from which we see that A,i = is one eigenvalue of Hq . 
In the single atom case the dressed state corre- 
sponding to the zero eigenvalue was of very special signif- 
icance as corresponding to a trapping (dark) state com- 
pletely decoupled from the fields |0, 0, [s^- Let us 
investigate this possibility in the multiatom case. 
If we represent the eigenvector |-D„) n by the column vec- 
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tor 



C2 



\D„ 



(27) 



we find by substituting into Eq. H24() that for A„ = the 
coefficients Cn for even n are all zero, whereas for odd n 
the coefficients are given by the recurrence relation 



C2fe+1 — 




1,2, 3,... (28) 



and ci is found from the normalization condition. 
The dressed state corresponding to the eigenvalue A„ 
can thus be written as 



\D)n = \D„)n 

= {cos9fY.(^^ )"(-tan0)^-|2^-), (29) 



fc=0 



where |2'') = |3), and 



tanO = 



O3 

^2 



(30) 



The collective dressed state (|29|l is a linear combination 
of the state |3) and of the states |2*''). The important 
feature of the state is that it does not contain the ex- 
cited states of the atoms and thus does not radiate. The 
dressed state is a stationary state of the Hamiltonian 
Hq describing the atoms driven by two coherent fields. 
Therefore, if nothing else is allowed to interact with this 
system, the state (|29|) will never change in time. 
We can write the multi-atom dressed state (I29|l in the 
basis of the single-atom dressed states Surprisely, 
we find that the state is of the form 



\D] 



N 



1*1) 



l*l> 



1*1; 



N, 



(31) 



which is a product of the single-atom trapping states 
|\E'i)j. Obviously, the state (|31|l is not entangled, which 
shows that trapping of the population in all of the atomic 
ground states is equally effective in destroying collective 
(entangled) properties of the system. Thus, the improve- 
ment of the CPT in the collective multiatom system, seen 
in Figs.01and[Sl does not arise from collective excitations 
of the dark state \D)i\f. 

We note in passing that the state in Eq. (|29|l is similar 
in form to that found by Mewes and Fleischhauer [l8l |. 
see also 0, > '^lio considered collective quantum mem- 
ories in three-level atoms driven by a classical field and 
a single-mode quantum field. The results of their work 
demonstrate that the dark states of the multiatom sys- 
tem are highly entangled states. However, the state \D)n 



which is the analog of the dark states found in is 
not entangled. Thus, a question arises: Why does the 
dark state \D)]s[ is not entangled? The reason is that 
the state \D)n is a linear superposition of all the collec- 
tive ground states, whereas the dark states considered in 
Ref. [I8j | is restricted to having involved a small number 
of the ground states corresponding to a small number of 
excitations fc <C iV. It is easily verified that if we limit 
the number of the states involved in the superposition 
(|29|l to k < N, then the resulting state cannot be writ- 
ten as a product of the single atomic states. Clearly, 
entanglement properties of the dark state \D)]\[ depend 
on the number of atoms involved in the interaction with 
the laser fields, that only for k < N the interaction can 
produce a dark state which is an entangled state. 
To find the explanation why the CPT in the collec- 
tive system interacting with the thermal field decoheres 
slower than the system of independent atoms, we intro- 
duce the interaction between the collective dressed states 
and the thermal field. This interaction leads to a distri- 
bution of the population, initially trapped in the dark 
state |-D)Ar, among the collective dressed states. Let us 
look at the evolution of the population of the state \D)n. 
Using the master equation H1U|) . we obtain the following 
equation of motion for the population of the state \D)]\[: 



PDD = ~'iNT2PDD + 2NTi {pD2 + PD3) 



(32) 



where puu is the population of the state \D)n and pD2 
are pns are populations of the following superposition 

states 

(1) ,T,(2) ,T,0") ,T,W 



2 /, 



i=i 

N 



i=i 

which differ in energy from the state \D)n by +fl and 
~rt, respectively. 

The states \D2)n and \D3)n are linear superpositions 
of the product states in which — 1 atoms are in state 
\'^i)j and one atom is in the state |^'2)j and |^'3)j, re- 
spectively. It is interesting to note from Eq. (|32(l that 
both spontaneous emission and the thermal field couple 
the state \D)n to only those states which differ in the ex- 
citation by one. Moreover, the transition rates between 
these states are N times larger than that for single atoms. 
This shows that the system is superradiant despite the 
fact that the state \D)pf is the product state of the single- 
atom trapping states. In addition, the collective decay 
rate of the radiators entering the state \D)n is larger 
than that describing the atoms escaping from it. Thus, 
the collective properties of the system are preserved due 
to the presence of the superposition states involving the 
single-atom states |\E'2)j and l^'s)^. 

In fact, the master equation H1U|I leads to a set of 
{N + 1){N + 2)/2 coupled equations of motion for the 
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■H Q 



2 3 
n 



FIG. 7: The ratio p^D2l pD2 as a function of n2 = na = n 
for different numbers of atoms: N = 2 (solid line), A'^ = 4 
(dashed line), A*' = 20 (short dashed line). 



populations of the collective dressed states. Fortunately, 
however, an explanation of the enhancement of the CPT 
effect, seen in Figs. 0] and |5| does not require a com- 
plete solution for the populations of the dressed states, 
ft is enough to consider only the population of the state 
\D2)n or IDS) AT. Thus, using Eqs. ® and (0, we can 
show that for = 1 the stationary population of the 
state \D2)m is simply equal to {R22), and for > 1 is 
given by the following expectation value 



PD2 



_ 1 

m 



{{R12R21 + R22 — Rii)Rii{Rii — 1) 
2)...(i?n-A^ + 2)), (34) 



which can be evaluated using the steady-state solution 

m- 

We also calculate the population of the state \D2)n in 
the case of independent atoms and find 



PD2 



r2 (Fi) 



N-l 



(Fi + 2F2) 



N ■ 



(35) 



To see the difference between the populations (|34ll and 
we study the ratio Pi)2/pD2- Figure [7| shows the ra- 
tio for different numbers of atoms. For iV = 1, the ratio 
is equal to one, but for N > 1 the ratio is smaller than 
one and decreases with N. This shows that the popula- 
tion of the collective states is larger than the population 
of the equivalent states of independent atoms. In other 
words, we may say that the capacity of the collective 
states is larger than the capacity of the equivalent states 
of independent atoms. 



The above analysis give us a simple physical interpreta- 
tion of the collective trapping effect. We may conclude 
that the improvement of trapping effect by multiatom 
system is due simply to the increased storage capacity 
of the collective (entangled) states compared with the 
storage capacity of the equivalent states of independent 
atoms. 

V. SUMMARY 

We have investigated the coherent population trapping 
effect in a collective system of three level atoms driven 
by two coherent laser fields and simultaneously coupled 
to the reservoir of a non-zero temperature. The thermal 
reservoir causes thermal decoherence which affects the 
trapping effect. We have shown that in a single atom 
there is no perfect CPT when both atomic transitions 
are affected by thermal decoherence. The perfect CPT 
may occur when only one of the two atomic transitions is 
affected by thermal decoherence. Extending the analysis 
to multi-atom systems, we have shown that the destruc- 
tive effect of the thermal decoherence on the CPT can 
be circumvented by the collective behavior of the atoms. 
Unlike the case of noninteracting atoms in which decoher- 
ence processes are independent of the number of atoms, 
we have found that the collective behavior of the atoms 
can substantially improve the trapping effect destroyed 
by the thermal decoherence. In the collective atomic sys- 
tem the trapping effect increases with increasing number 
of atoms. If number of atoms is large enough, an almost 
complete CPT is observed even at high temperatures of 
the reservoir. This feature is explained in terms of the 
semiclassical dressed atom model. We have shown that 
the improvement of the CPT trapping in the multiatom 
system arises from the presence of collective (entangled) 
states whose capacity of storage of the atomic population 
is larger than the corresponding states of independent 
atoms. 
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